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SUMMARY 

A study was made of the .lam.1.nar conrpressi'ble Taoundary layer induced 
by a shock wave advancing into a stationary fluid bounded by a wall. For 
weak shock waves ^ the boimdary layer is identical with that which occurs 
when an infinite wall is impid-sively set into uniform, motion (Raylei^ 
problem) . A numerical solution was required for strong shocks. 

Velocity and temperature profiles, recovery factors, and skin- 
friction and heat-transfer coefficients are tabulated for a wide range of 
shock strengths. 


INTRODUCTION 

If a shock wave advances into a stationary fluid bounded by a wall, 
a boundary-layer flow is established along the wall behind the shock. 

This boundary layer is often important in studies of phenomena involving 
nonstationary shock waves. In a shock tube, for example, this boundary 
layer acts to attenuate the strength of the shock which propagates throu^ 
the low-pressure side of the tube (refs. 1 and 2) . If the shock tube is 
used as an aerodynamic wind tunnel, the test time available may depend, 
for long shock tubes, on the time it takes the boundary layer to Introduce 
nonuniformities in the test section. 

Another exanple of a shock-generated Doundary layer occurs when a 
combustible mlxtvire is ignited within a tube. In this case, a shock wave, 
followed by a flame front, is observed, as discussed in references 3 to 
5. The shock wave is particularly strong when ignition occiirs at a closed 
end. For long tubes, the progress of the flame front will be related to 
tlie boundary-layer development behind the shock. Since flame speed is in- 
creased by fluid txirbulence, a transition from laminar to turbulent flow 
will accelerate the flame. T3ius, the boundary layer may play a role in 
the acceleration of a low-speed flame to a detonation wave in a long tube 
(ref . 4) . 
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!Hie boundary layer ’behind the shock waa studied in references 1 and 
2. Both, papers were primarily concerned with shock -wave attenuation. 

In reference 1, an approodLmate solution for the boundary layer was esti- 
mated 'by referri.ng to the flow induced within a clrcrular cylinder that 
is Impulsively set into uniform translation. The validity of this anal- 
ogy was not established. In reference 2, the correct boundary-layer equa- 
tions were considered. These were integrated with a REAC (Beeves Elec- 
tronic Analog Computer) . "Values of skin-friction and heat-transfer coef- 
ficients were presented. However, no velocity or temperature profiles 
were reported. Because of the growing interest in phenomena related to 
these sbock-lnduced boundary layers, it was felt that a more detailed 
and more accurate study of this boundary- layer problem was warranted. 

Such a study was conducted at the MCA Lewis laboratory and the results 

are presented herein.^ 

In the following sections, the laminar compressible boundary layer 
behind a shock wave advancing into a stationary fluid, bounded by a wall, 
is analyzed. For weak shocks, an analytical perturbation solution is pre- 
sented. Numerical resvilts for velocity and temperature profiles and heat- 
transfer and skin-friction coefficients are tabulated, covering the range 
from weak to strong shocks. The numerical results are correct to four 
decimal places. 


AMLYSIS 

Coordinate ^sterns 

A shock wave of constant strength, is considered to move, parallel 
to a wall, into a stationary fluid. Let (x,^ be a coordinate system 
fixed in_respect_to the wall and let u ’and ^ be velocities parallel 
to the X and y coordinates, respectively, as indicated in figure l(a) . 
The flow is unsteady in this coordinate system. Let (x)'y) represent a 
coordinate system moving with the shockwave (fig. 1(b)). The velocities 
parallel to the x and y coordinates are denoted "by u and v. In 
this coordinate system, the flow is steady. 

__Assume that, at time t := 0, the two coordinate systems coincide. _ 

If u is the velocity of the shock wave relative to the wall, then x 
s 

and X are related by x = x - u^t. Similarly, the axial velocities are 
related by u = u - u^. Note that the wall moves with velocity u^ = -u^ 
in the steady coordinate system. 


^TEhe writer's present Interest 
problem was stimulated ’by a private 
Cornell University, who is studying 
with shock tubes. 


in the shock-induced boundary- layer 
communication from Prof. H. Rott of 
heat-transfer problems associated 
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Bomdazy-Layer Equations 


The Prandtl houndBry-layer equations apply for the flow in the vicin- 
ity of the wall (except at the base of the shock, where the boundaiy-layer 
assumptions break down) . Assuming the flow to be laminar and dp/dx = 0, 
the equations of motion are, for x > 0, 


^ = 0 (continuity) 




u + V 
ox 


pc. 




3 ^ 


+ V 






p = pRT (state) 


J 


CD 


The additional symbols are defined in appendix A. The boundary conditions 
for X > 0 are 


U(X,0) = -Ug 

u(x,-) 

v(x,0) == 0 

T(x,») 

o 

II 




= T. 




( 2 ) 


These are the usual bomdary conditions, except that the fluid at the wall 
moves with velocity u(x,0) s u^ = -u^ in order to satisfy the condition 

of zero slip at the wall. It will be assimied that_the wall temperature 
Tjj. is constant. The magnitudes of Ug, Tg, and Ug depend on the shock 

strength and can be found from the normal shock relations quoted in ap- 
pendix B. The ratio ii^/Ug Increases from a value of 1 for a very weak 

shock to a value of (r + l)/(r - l) for a very strong shock wave. Thus, 
in the steady coordinate system, the u velocities in the boundary layer 
have a maximum at the wall and decrease monotonically to the value in the 
free stream (as indicated in fig. 1(b)). 


Transformation . - Equations (l) and (2) are transformed to a system 
of ordinary differential equations with the methods of references 6 and 7. 


From the continuity equation, a stream function i|r exists such that 
3t/3y = pu/p^, = pv/py.. Following reference 6, a similarity parame- 

ter Tj is defined according to the relation 
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n = 




2 XV, 


V 



and "the stream function is written as 




(3) 


( 4 ) 


Note that f’ = u/ug. As in reference 1, the viscosity is assumed to vary 
linearly with temperature . If the viscosity is referenced to the wall 
value, 



Substitution of equations (3), (4), and (S) into the momentum equa- 
tion yields 

f * • • + ff" = 0 (6) 


with the boundary conditions 


f(0) = 0 


f«(o) = 

f (-) = 1^ 


(7) 


Equation (6) is the familiar Blasliis differential equation. However, the 
tangential velocity boundary condition at the wall (f '(O) = u^/u^) is dif- 
ferent from the aero value usually encountered in studies of viscous flow- 
past a semi -infinite plate. Numerical integration is required except for 
the limiting case of a weak shock, [(u^/u^) - 1] «1, for which an 

analytical perturbation solution is possible. 


For T a function of ti only, the energy equation becomes 

T" + cjfT» = -a(r - 1) M|(f")^ C8) 

assuming that -the Erandtl number a is constant. Since equation (8) is 
linear, the general solution for T can be expressed as -the linear super- 
position of the solution for zero heat transfer plus the effect cd* heat 
transfer. That is, T/Tq can be expressed in -the form 
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where r(Ti) satisfies 


r" + cr fr* 


r(-) 



( 10 ) 


and s(t\) satisfies 


s" + ff fs* =0 
s(0) = Ij s(») = 0 


( 11 ) 


Note that "ll^-i^glAe 1® Mach, number of the external 

flow relative to the wall. For an insxxlated wall, the coefficient of 
s (in eq.. (9)) equals zero, so that the wall temperature is 





- 2 

Mg r(0) 


( 12 ) 


Thus, r(0) is a recovery factor based on the Mach number of the external 

flow relative to the wall. Equation (9), in terms of I . , is 

w,l 








Equations (lO) and (ll) can be expressed in quadratiire form 


(12a) 


r 




[f ”(€)]'" d€ 



[f"(e)3^"'^ d9 


(13) 
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s = « 

For O' =• these equations can he integrated to yield 



e 


The solutions for other Prandtl nianhers are discussed later. 


(14) 


(15) 


(16) 


An alternate system of equations is described in appendix C. 

Relation between y and ti . - For x constant, the relation be- 
tween y and q is, from equation (3), 


y 






dn 


Substitution of equation (l2a) into equation (17) yields 


(17) 



(18) 


For cr = 1, equations (l5) and (16) can be substituted into equation (l8) 
with -the following res-ult: 
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y A o 


1 


2 % 



(f - n) + 


IL^m2 

2 e 


u 


^ (f - n) + f(i - f •) + f"(o) - f"| > (19) 


If the wall is insulated^ equation. (l9) hecomes 


'1 e 




(f - - 1^ + f(l - f ') + f"(0) - 


f" 


y/\ 2 


xv_ 


w 


1 + 


r - 1 




w 




( 20 ) 

Equations (l8) to (20) axe useful for obtaining velocity and temperature 
profiles in terms of y rather than t\. 


Perturbation Solution for 



«1 


If the shock wave is weak, (u^/ug) - 1 is small, and a perturbation 
analysis in terms of this parameter is possible. Let 


f 



(21a) 


r 




s 



(21b) 


(21c) 


Substituting equation (21a) into equations fs^ and (7) and equating coef- 
ficients of [(Uj^/u^) - 1]® [(u^,/Uq) - 1]^ yield, respectively. 
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^o'*' W =0 
foCo) = ^6(0) = 1 

^O^l" ■*■ ^l^o" ^ 0 
f^(0) = 0) f^(0) = 1; f^(-) = 0_ 
Integration of equations (22) and (23) gives 

fo = ^ 


( 22 ) 


(23) 




\ 

n } 

J 


(24) 


so that 


f = Ti + /Ti erfcf 


f . . 1 ^ [erf=(-^J(^ - ij * o(^ - ij 

^■' - -V! - °fe - ^ 

Substituting the preceding solution for f into equations (lO) 

.0 


oN 


^(25) 

J 

and (11) 


and equating coefficients of 




yield 




ro + anr^ = 0 e 

^6^^) = ^o(*) = ° 

Sq" + cniS(^ = 0 

SqCo) = Sq^*^ " ° 


(26) 


( 27 ) 
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Integration, of equation (26) yields (from results of ref. 8) 


sa! as + ofiis - 1 


4/ C \V2 

-*( rvr ) “^(-2 

J 


2 sin^e 


u. 


(28) 


The Integral of equation (27) is 


s = erfc 




(29) 


From equations (28) and (29) , 


“■( 0 ) = I (2^) ^ (1 - Vi) + °(s ■ ^ 




(31) 


where r(0) is the recovery factor and s’(o) is used to calculate heat 
transfer. Equations (22) to (31) (neglecting the higher-order terms) 
wply if an infinite plate is started, impulsively, with velocity 
[(ujf/ug) - 1 ]. The latter, often termed the "Eaylei^ problem, " has been 

much discussed. The zero-order solution indicated by equations (28) and 
( 30 ) was obtained in reference 8 in a study of the aerodynamic heating 
associated with the Rayleigh problem. 

In reference 1, the boundary layer behind a shock wave advancing into 
a stationary fluid was estimated by analogy with the Rayleigh problem. 

The work of the present section shows that this approach is exact for weak 
shocks . 


Numerical Solution 

For other than weak shock waves, a niJtmerical solution of equations 
(6), ( 10 ), and (ll) is required. Th3 s was obtained by lynn U. Albers with 
the use of an IBM card-programmed electronic calculator. The integration 
technique is described in appendix B of reference 9. The results, correct 
to foiir decimal places, are tabulated in table I. Values of f, f', and 
f" are given for u^/ug of 1.5, 2, 3, 4, 5, and 6, while values of r. 
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rS s, s', r dti, and s dij are presented for of 2 , 4, 

and 6 (vrlth. cr = 0.72). The data of talle I pl\is the perturbation solu- 
tion of the previous section define solutions covering the range from very 
weak shocks to very strong shocks. 


EESULTS AND DISCUSSION 


In the steady coordinate system, the boundary- layer similarity param- 


eter is 


1 ^ 


Tw 


!Eransformlng to the unsteady coordinate 


system according to the relation x = x + u^^t and considering station 


X = 0 gives Ti = 



dy. Thus, the boundary layer 


behind the shock wave has features of a Blasius type flow or a 'Rayl ej gh 
type flow, depending on whether the observer is stationary with respect 
to the shock wave, or wall, respectively. 


The Raylel^ viewpoint is used herein to correlate and discuss the 
ntanerical data. That is, attention is fixed at station x » 0, and the 
boundary-layer development for t > 0 is considered. The velocity which 
characterizes the boundary-layer development is u -u^. aimllarly, a 
characteristic length is (u^. - Ug) t, which is the distance a particle in 
the free stream moves relative to the wall in time t. The form of the 
Reynolds number used herein can then be defined as 


Re = 


V 

w 


(32) 


Boundary-Layer Profiles 

The parameter t(u^/ug) - f ']/[(xi^^/ug) - ij varies from a value of 
zero at the wall to a value of 1 at the edge of the boundary layer. A 
boundary-layer thickness 6 may be defined as the value of y correspond- 
ing to [(u^^/ug) - f']/[(iv/^e) ■ =0.99. Values of the boundary-layer- 

thickness parameter have been computed assuming y = 1.4 and 

an insulated plate, for cr of 0.72 and 1.0 (using eqs. (l8) and (20)). 
These are tabulated in table H. It is seen that (B//i/'^)i increases 

with increasing u^/ug . As escpected , the values for n = 1 are greater 
than those for cr » 0.72 (for u^/ug ^ l) . This is due to the fact 
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that the larger recovery factor of the former leads to hi^r temperatures 
near the wall and therefoire a larger value of y corresponding to a 
given Ti . 

Velocity profiles for of 1 (limiting case of very weak 

shocks) and 6 are plotted in figure 2 . Curves for intermediate values 
of Uy/Ug lie smoothly between the curves in the figure. No marked 

departure from an error function type velocity profile is indicated. 


Skin Friction and Heat Transfer 


The shear stress at the wall is given hy 

Because of the coordinate system used, is negative. If (u^ - Ug) is 
used as a reference velocity, a positive local friction coefficient can 

he defined as c^ = - Ug)^. Then, using the Beynolds number 

as defined hy equation (32) , 


= 


■A[Zt"iO) 

_ A 


(33) 


Values of are tabulated in table II . These vary from the Ray- 
lei^ value of 1.128 at ^ 0.935 at h^/^g = 6. The corre- 
sponding value of Cj/\/Re for incompressible flow past a semi -infinite 
plate (Blasius problem) is 0.664. 

The heat transferred into the fluid from a unit area of the wall, 
per unit time, is 


1 = = A V? ^ 

Defining a heat -transfer coefficient h by h = (T^ - s- 

Nusselt number by the relation Nu = h(u^ - Ug) t/k^ permits the Kusselt 
number to be written as 
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The relation between skin friction and heat transfer can be exjiressed in 
terms of a Reynolds analogy parameter c^Re/Su. The factor VcT ls_. 

included, since, for u^j./ug = 1, /{a c^e/Nu is a constant. For a = 1, 
the parameter equals 2 at all values of while for a = 0.72, the 

parameter increases from 2 at =1 to 2.07 at ~ 

The recovery factor for u^/ug =1 1 b given by equation (30) . Evalu- 
ating this equation for a = 0.72 gives r(o) = 0.885. The recovery 
factor for a = 0.72 increases with increasing to a value of 0.920 

at = 6. These compare with the value 0.845 for flow past a seml- 

infinlte plate at 0 = 0.72. 

Thus, for the range of investigated, the numerical results 

for skin friction, heat transfer, and recovery factor (in terms of the 
parameters defined herein) depart relatively less from Raylel^ (u^/ug * l) 

values than from the Blaslus values for equivalent flows past a seml- 
inf inite flat plate . 


CONCLUDIRG REMARKS 

The laminar boundary layer behind a shock wave advancing into a sta- 
tionary flviid, bounded by a wall, has been determined. Various boundary- 
layer parameters have been tabulated for several shock strengths. 

With increasing Reynolds numbers, the laminar boxindary layer behind 
the shock will become unstable, and transition to turbulent flow will 
ultimately occur. A theoretical study of the stability of this laminar 
boundary layer would be of interest. Shock-tube experiments might provide 
criteria defining the transition to turbulent boundary- layer flow as well 
as the characteristics of the turbxilent boundary layer. At present, little 
is known about the structure of such turbulent boundary layers, and it is 
felt that some experimental data should be available before an analytical 
study is attempted. 


Lewis Fli^t Propulsion Laboratory 

National Advisory Committee for Aeronautics 
Cleveland, Ohio, December 10, 1954 
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APPENDIX A 


SYMBOLS 


The following symbols are used in this report: 
a speed of sound 

p (u 




local sMn-frlction coefficient. 




erf(x) 

erfc(x) 

f 

h 

k 

M 

Me 

Nu 

P 

q. 

R 

Re 

r 

r(0) 

s 

T 

t 

u,v 


specific heat at constant pressure 

error function, (2/ e~^ dy 

complementary error function, 1 - erf(x) 
function of r\ defined hy eq. (4) 
heat-transfer coefficient, (T^. - T^^j^)/q 

thermal conductivity 
Mach number 

^e/^e 

Nusselt numher, h(u^ - Ug)t/l^ 
pressure 

local rate of heat transfer 
gas constant 

Reynolds number, (u^ - Ug)^t/v^ 
function defined by eq. (lO) 
recovery factor 
function defined by eq. (ll) 
static temperature (abs) 
time 

velocities parallel to x and y coordinates, respectively 
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velocities parallel to x and y coordinates, respectively 
Ug velocity of shock vave relative to wall 

x,y coordinates parallel to and normal to wall, respectively, and 

moving with shock’’wave (fig. 1(h)) 

x,y coordinates parallel to and normal to wall, respectively, and 

stationary with respect to wall (fig. 1(a)) 

X ratio of specific heats 

6 value of y for which (u^^Ug - u/ug)/(u^/ug - l) »* 0.99 

(i.e., boundary- layer thickness) 

T) variable defined by eq.. (3) 

ilg value of T] at y ■ 6 

|i coefficient of viscosity 

V kinematic viscosity 

p mass density 

cr Prandtl number, pCp/k 

local shearing stress at wall 

t stream function 

Subscripts: 

b undisturbed flow in front of shock wave (appendix B) 

e flow external to boundary layer 

i insulated-wall problem 

conditions at wall (y » y = O) 


w 
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APPENDIX B 
SHOCK RELATIOI^ 

Eor convenience, scane useful shock-wave relations are noted herein: 


U^, EU^ ^ 

y 

®b» 

f ^e/ He 




Consider flow in a steady coordinate system. Let subscript h designate 
undisturbed flow in front of the shock wave; and subscript e, the flow 
behind the shock wave and external to the boundary layer. Note that 
Ujj e u^j so that Mjj = %/% ” '^/®b* 








(r-i)i4 + 2 

6m2 


+ 5 


for T at 1.4 




(Bl) 








(r+1) ^ - (r-i) 


6^-1 


for X * 1*^ 




(B2) 


T 


e 

Tb 


r N Sr 

III 

e 

- (r-1) 

Ug 

|ir+i) - 

(r-x) §] 

6 

ue " 

for r 





6 

'V ' 





1.4 




(B3) 
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APPENDIX C 


ALTEREAIEE FORMULATION OF BOUNDARY -LAYER PROBLEM 

The transformation based on equation (s) leads to a system of equa- 
tions identical to that which arises in boundary- layer stupes of the 
flow past a semi -infinite flat plate, except for the non-zero velocity 
boundary condition at the wall. This sytem was convenient for numerical 
computation, since a variety of flat-plate boundary-layer problems had 
previously been solved at the NACA Lewis laboratory, and ■Ube card- 
programming for the IBM electronic calculator was already established. 

An alternate system can be obtained by normalizing the momentum- 
equation boundary conditions. That is, the parameter appears in 

the differential equation rather than in the boundary conditions. Such 
a system leads more dlfectly to the correct physical parameters of the 
problem and is described by the following equations. 

Define ^ «nd g(^) according to the relations 



Using equations (Cl) and (C2) , equations (6), (7), (lO), and (ll) become 




G 
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a 







ds 


0 




s(0) *S Ij s(i») s 0 




(C5) 


Hote that, for x = 0, ^ = l/y2v^t, -which is the correct parameter for 
the Rayleigh problem. Also dg/d^ = [(u^/ug) - f *]y/[Cuy/ug) - 1], which 

is a normalized velocity. Finally, the reduction of the system to a 
Rayleigh problem, for Uy./ug * 1, is more apparent from equations (C3) 

to (C5) than from equations (6) to (ll) . 
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KABia I. - LU02JAR BOUBMHY UUCKH HBUBD SHOW WAVS 
(a.) Solution of voaontua (quntlon for > 1,3, 3.0 and 5,0 


H 

- 1.5 

V“8 ” =‘° 

'V"e ' s.o 

■■ 


f 

-f" 

f 

f ' 

-f’ 

f 


-f" 


^WSRSffi^l 

Bl TlSIi^l 

0.4 5 7 0 

OO 00 0 

3,0 0 00 

a3 973 

043000' 

&0 0 00 

397 3 6 

.1 


kISkH 

.4 544 

■88 8 0 

a? 614 

a3620 

.4 70 3 

4.4 076 

5,83 0 9 

.2 

3909 

1.4 09 3 

.4 445 

.5525 

34 3 96 

a2 64 9 

.8924 

38 419 

344 7 0 

-.3 


1.3 656 

.4 288 

.794 4 

3.3 1 0 1 

ai 170 

1.8 4 0 3 

33 2 37 

4.89 64 

.4 


1,3 3 38 

.4 080 

1.0 15 1 

2J. 0 74 

L9336 

LS492 

aa 6 56 

4.2 5 7 3 

.5 


1,2 84 2 

.3 8 31 

1.2 16 6 

1.9 241 

L7 39 1 

1.8 1 5 5 

a4 729 

35968 

.6 


1J3 47 2 

■3651 

1.4 00 6 

1.7 618 

LSI 68 

ao 4 59 

ai 4 62 

Z.9 64 4 

•7 


1.2 1 3 3 

.3 251 

1.56 9 6 

1,6 807 

L3073 

2.2466 

L8 78 0 

&3 913 

■a 

1.0 6 37 

1382 3 

.29 41 

1.7 2 5 4 

Lao 00 

L10 87 

2.4 2 33 

L6 645 

L8 93 0 

.9 

1J.0O6 

14. 54 4 

.2 6 28 

1.0 7 0 a 

1.3 9 84 

.9 26 1 

2.5 B ). 0 

L4 968 

L4 73 7 

1.0 

1,2 94 7 

1.1 2 9 7 

.2 383 

3.0 0 5 7 

1.3141 

,7 629 

3.7 2 3 9 

L3673 

L130 3 

1.1 

1,4066 

1407 9 

.2 0 29 

2.1 3 3 S 

1.3 4 51 

.6 303 

a8555 

L3 68 5 

,8551 

1.2 

1.5 164 

108 90 

.176 3 

ae 55 1 

JL18 94 

,4 98 0 

2.9784 

L194 3 

.63 8 7 

1.3 

1,6 24 5 

1.0 7 2 8 

4 496 

a3 71 a 

1.14 49 

.3 95 3 

3.09 50 

LI 391 

.4714 

L4 

1,7 3 1 0 

1.0 5 9 0 

.136 7 

^4 84 4 

1.1 0 97 

.3100 

3,3067 

L0 98 6 

.344 0 

1.5 
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TABLE II. - SKIN-FRICTION AND HEAT-TRANSFER COEFFICIENTS 
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(b) Cciordinate system fixed with respect to shock wave. 

Figure 1. - Coordinate systems used to study boundary layer behind shock advancing into 
stationary fluid. 
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Figure 2. - Boundary-layer velocity profile, 
fers to limiting case of very weak shocks . 
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